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ABSTRACT 

How thick should the sample be for a transmission THz-TDS measurement? Should the sample be as thick 
as possible? The answer is 'no'. Although more thickness allows T-rays to interact more with bulk material, 
SNR rolls off with thickness due to signal attenuation. Then, should the sample be extremely thin? Again, the 
answer is 'no'. A sample that is too thin renders itself nearly invisible to T-rays, in such a way that the system 
can hardly sense the difference between the sample and a free space path. So, where is the optimal boundary 
between 'too thick' and 'too thin' ? The trade-off is analysed and revealed in this paper, where our approach is 
to find the optimal thickness that results in the minimal variance of measured optical constants. 

Keywords: THz-TDS, T-rays, terahertz, optimum thickness 

1. INTRODUCTION 

Terahcrz time-domain spectroscopy (THz-TDS) is a promising tool in determining the complex response of 
materials in the T-ray frequency regime, bounded between 0.1 and 10 THz. The system generates and detects 
coherent broadband pulses, usually on the basis of ultrashort laser excitation. Recent developments of the system 
mainly aims for wider bandwidth 1 ' 2 and higher SNR, 3 but not the uncertainty of the signal. It is known that in 
addition to the noise in optics 4 and electronics, 5 the mechanical drift 6 results in fluctuation of measured signals 
over time. 

Not many methods arc available in reduction of the signal uncertainty. Increasing the signal strength and 
using optical chopper can lessen the contribution of noise from electronic parts, but cannot solve the problem 
from mechanical drift and optical noise. Repeating the same measurement a number of times, in practice, allow 
the effect from mechanical drift to manifest in measured signals. Another different approach in reduction of the 
uncertainty in signals is hidden in the sample under measurement. 

It is well known to every THz-TDS experimentalist that, in the transmission-mode spectroscopy, too thick of 
a sample with considerable bulk absorption can significantly reduce the signal power and increase the uncertainty 
of measurement. In response to that awareness, a sample is usually made very thin, while its strength and shape 
are well retained. However, too thin of a sample can also cause problem, as the system might not be sufficiently 
sensitive to detect the change in amplitude and phase of the signal. In this paper, the centreline between the two 
extremes is proposed. The determination of this line is on the basis of minimum uncertainty in measurement. 

The article is organised as follows: In Section 2, an analytical model relating the variance in signals to the 
variance in optical constants is introduced. This model leads to the optimisation of the sample thickness by 
minimisation of the variance in optical constants, as shown in Section 3. An analytical formula for optimum 
thickness, as the outcome of the optimisation, is verified by the experiments with various materials in Section 4. 
The usage of the formula is discussed in Section 5. The article ends with conclusion in Section 6. 
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2. VARIANCE IN OPTICAL CONSTANTS INFLUENCED BY NOISE 



The T-ray amplitude is prone to variation induced by many sources of random and systematic errors. The 
sources of random error include laser intensity fluctuation, optical and electronic noise, jitter in the delay stage, 
etc., whereas the sources of systematic error include registration noise, mechanical drift, etc. The variation in the 
amplitude may embrace the effects from inhomogencity in a sample or among samples, if the sample is displaced 
or replaced with nominally identical samples during several measurements. Considered here is the amplitude 
variance model, which unites all these errors and assumes a normal probability distribution. The influence of 
this amplitude variance on the extracted optical constants is shown in this Section. 

Given that the sample under measurement has parallel and polished surfaces, and the angle of incidence of 
the incoming T-ray beam is normal to the surfaces, the transmissive transfer function of the sample is expressed 
as 
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where -E re f(w) and E sam (oj) is the time domain reference and sample signals, I is the sample thickness, n(ui) and 
k(uj) are the refractive index and the extinction coefficient of the sample, Uq is the refractive index of air, and r 
and t' are the transmission coefficients at the sample interfaces. 

Influenced by the variance in measured signals, the variances in the optical constants, derived based on 
Equation 1, read, 
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Here, s 2 E ^ { {k) and s E ^(k) are the variances associated with the reference and sample signals, respectively; k is 
the temporal index and r is the sampling interval, and thus kr is the time. The summation is carried out over 
the time duration of a recorded T-ray signal. In the equations, all the parameters are calculated at their mean 
value. The proposed model in Equation 2 is successfully validated with Monte Carlo method. 



3. OPTIMISATION OF THE SAMPLE THICKNESS 

From Equations 2a and 2b in Section 2, it can be inferred that four major variables, n(ui), k(u>), I, and |.B re f|, 
govern the amplitude-related variance, E or s 2 K E . Optimising one of these parameters might help reduce the 
variance. Two of them, n{uo) and k(oj), are intrinsic to materials, and thus cannot be optimised. The signal 
strength, |-E re f|, with no doubt, must be as large as possible providing no damage to the sample to minimise 
the variance. The sample thickness, /, is a flexible factor, which may lend itself to optimisation. 

The intricate relation between the sample thickness and the standard deviation in the refractive index, s n E , is 
simulated and demonstrated via a contour plot in Figure 1. Apparently, at every frequency, there is an optimum 



sample thickness that gives the lowest s n ,£. Figure 2 reveals the magnitude of s n ,E and s k ,e as a function of the 
thickness, estimated at three different frequencies. The optimum thicknesses for the simulated sample at these 
frequencies are grouped around 300 /im to 1 mm. Moving towards a thicker sample by an order of magnitude 
sees an increment of the standard deviation by three orders. Likewise, moving towards a thinner sample by an 
order of magnitude results in one order stronger standard deviation. Selecting the sample thickness to match its 
optimality is therefore advantageous. 

According to the transfer function in Equation 1, 
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Substituting Equation 4 into Equation 2a and reassigning the notation gives 



where 



C = zm^w and = (6) 

Now we are able to minimise E (u>) with respect to the thickness I. By taking the derivative of Equation 5 
with respect to I, we arrive at 
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Figure 1. Standard deviation in the refractive index against the sample thickness and frequency. The contours represent 
the simulated magnitude of s n ,E, which is in fact comparable to that of s Kt E. The optical constants are set ton - jn — 
3.0 — O.lj (k — 0.1 equals a — 41.9 cm -1 at 1 THz). An additive white Gaussian noise limits the maximum dynamic 
range of the reference spectrum to 40 dB. The profile of the optimum thickness determined numerically appears jagged 
due to the discrete nature of the simulation. 
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Figure 2. Standard deviation in the optical constants against the sample thickness. The plots are cross-sectional profiles 
of Figure 1, picked up at 0.5, 1.0, and 1.5 THz. The optimum thicknesses, where the standard deviation is minimal, are 
indicated by arrowheads. By moving towards the thicker sample, the standard deviation rapidly increases to the point 
comparable to the optical constants' values. The standard deviation at high frequency is more sensitive to the thickness 
increment, as the T-ray magnitude at high frequency is relatively low. 



Using the above equation, the amplitude variances in the reference and sample spectra can be related via the 
transfer function, 
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Substituting Equation 8 into Equation 7 and equating to zero gives, 



4 ^V" 2 ^) ^rcfM - 2^U rcf ( W ) = 0, (9) 

and further manipulation yields the optimum thickness: 
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Optimisation of the sample thickness, /, by starting from Equation 2b also delivers the same outcome. Notice 
that the optimum thickness, Z op t, relies on neither the index of refraction, n{u), nor the signal magnitude, 

In Table 1, some common dielectric materials are determined for their optimal thickness at 0.5, 1.0, and 
1.5 THz, according to Equation 10. It appears that the resultant thickness calculated from this equation is also 
agreeable with the result found numerically in, 10 in which the optimisation is for the thickness difference in 
dual-thickness THz-TDS. For example, at 1 THz, the optimum thickness for dioxane, which has k — 0.013 or 
a = 5.45 cm , is reported to be 4 mm, and the optimum thickness for water, k — 0.478 or a = 200 cm , is 
100 /im. Using the same parameters, the derived analytical expression here estimates the thicknesses for dioxane 
and water at 3.7 mm and 100 /im, respectively. 



Table 1. Optimum sample thickness for some common materials. The optimum thickness for THz-TDS measurement 
is determined using Equation 10. The absorption coefficients, measured at room temperature, are taken from various 
sources: water 7 ; PMMA, u.h.m.w PE, TPX 8 ; HDPE. 9 HDPE — u.h.m.w is ultra-high molecular weight polyethylene. 
Note that the absorption coefficients can vary largely from sample to sample, in particular for plastics. Empty entries 
mean that data are unavailable. 
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4. EXPERIMENTS AND RESULTS 

The measurements are carried out with a fibre-coupled T-ray system from Picomctrix. The system generates 
the pulsed T-ray radiation spanning 0.05 to 2 THz, with the maximum dynamic range of 30 dB. The reference 
measurements are made in two batches, ten scans each, before and after every single batch of the sample 
measurements. The similarity of the variances of the two batches confirms no significant change to hardware 
during sample measurement. 

4.1. Polyvinyl chloride: PVC 

The measurements are carried out with normal-grade PVC, which is preformed in a rod shape. The rod with a 
diameter of 50 mm is cut into four cylindrical segments, with the length of 1 (0.9672), 10 (9.9772), 20 (19.8694), 
and 50 mm (the values in parentheses are 5-point averaged). The surfaces of the segments arc well-treated to 
minimise the scattering effect. Each segment is measured with a collimated beam from a THz-TDS system in 
the axial direction ten scans, and each scan is made after the previous scan by 30 sec. 

4.1.1. Optical constants 

The optical constants of PVC are extracted from the measurement of the 10- mm sample, as it provides the 
widest bandwidth among the rest, and its thickness is more uniform than the 1-mm sample. The phase at low 
frequency band is extrapolated from the phase value between 0.05 and 0.1 THz. Figure 3 shows the extracted 
optical constants. Clearly, the refractive index is nearly constant at 1.63, consistent with the published value at 
1.67. 11 The absorption coefficient increases quadratically, also in accord with the published value. 11 

4.1.2. Optimum thickness 

The optimum thickness for PVC, which yields the lowest variance in the measured optical constants, is determined 
from the absorption coefficient using the proposed model in Equation 10. However, the absorption coefficient 
contains error from noise, water-vapour absorption, and Fabry-Perot fringes. By assuming that the absorption 
coefficient is well described by a parametric function, the measured coefficient is initially smoothed by a second- 
order polynomial, which also enables extrapolation of the coefficient towards a higher frequency range. Figure 4 
illustrates the optimum thickness determined directly from the measured coefficient and from the fitting model. 
It can be seen that at frequencies around 1.0 THz, the sample thickness of 1 mm would provide the lowest 
variance of the optical constants. In addition, at low frequencies the optimum thickness increases by around one 
order of magnitude. 
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Figure 3. Optical constants of PVC. The constants are calculated from the signal probing the 10-mm-thick sample. The 
available bandwidth of the measurement is from 0.05 to 0.55 THz. The fringes at low frequencies are from Fabry-Perot 
reflections. 




Figure 4. Optimum thickness for PVC. The optimum thickness is determined from the measured absorption coefficient. 
The solid line is the second-order polynomial fitting to the absorption coefficient between 0.05 and 0.55 THz. The straight 
dotted lines indicate the thicknesses of 1, 10, 20, and 50 mm, available from the measurements. 
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Figure 5. Standard deviations in the optical constants of PVC. Each standard deviation is determined from ten time- 
resolved signals. The scatter plots represent the raw profiles obtained from Equation 2, whereas the solid lines are from 
an analytical function, exp(aia; 6 + 02a; 5 + . . . + agx + ay), fitting to the scatter plots. The missing part of the profiles at 
a high frequency range corresponds to the low SNR portion of the measured spectra 

4.1.3. Standard deviation 

Figure 5 shows the standard deviations of the optical constants for the four PVC samples, in terms of the 
unprocessed scatter plots and the fitting parametric curves. It is evident that at 0.05-0.2 THz, the 50-mm-thick 
sample provides the lowest standard deviation among the four samples; at 0.2-0.25 THz, the 20-mm-thick sample; 
at 0.25-0.6 THz, the 10-mm-thick sample; and above 0.6 THz, the 1-mm-thick sample. This optimum relation 
is in perfect agreement with the prediction in Figure 4, which is derived from the proposed optimal-thickness 
model. An example of the improvement in measurement can be observed by comparing the standard deviations 
of the 1-mm-thick and 50-mm-thick sample. At around 0.1 THz the value for the thicker sample is «4xl0 -5 , 
and that for the thinner sample is «2xl0~ 3 , or the improvement of the standard deviation is by almost two 
orders of magnitude. 

4.2. Polyethylene: PE 

The PE pellets used in the experiment are prepared by pressing amounts of pure ultra-high molecular weight 
(u.h.m.w.) polyethylene powder in a pellet press. The obtainable pellets' surfaces are an optical grade, and thus 
the surface scattering is negligible. The available thicknesses of the pellets are as follows: 4.5 (4.542), 5.4 (5.414), 
8 (8.066), 9 (9.142) mm. Note that a thicker sample is not available due to the pellet's instability. Each of the 
PE pellets is measured with a focused THz beam ten scans, and each scan is delayed from the previous scan by 
30 sec. 

4.2.1. Optical constants 

The optical constants of PE shown in Figure 6 are extracted from the 9-mm sample, as it causes obvious change 
to the signal. In the extraction process, the phase spectrum at the low frequencies is extrapolation from the 
phase at 0.1 to 1 THz. The average index of refraction is at 1.465, and the absorption coefficient is below 1 cm -1 
up to 2 THz. The measured absorption coefficient is in agreement with the value published in, 8 whereas the 
published value of refractive index for u.h.m.w. PE is not available. 

4.2.2. Optimum thickness 

Shown in Figure 7, the optimum thickness of PE is obtained by using Equation 10. As PE is virtually non- 
absorptive at low frequencies, no optimum thickness below 0.4 THz is determinable. In addition, the exceptional 
transparency of PE causes considerably large optimum thickness, i.e. higher than 10 mm at all frequencies. Thus, 
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Figure 6. Optical constants of u.h.m.w. PE. The constants are calculated from the signal probing the 9-mm-thick sample. 
The available bandwidth of the measurement is from 0.05 to 1.8 THz. The water-vapour resonances are observed at 1.1, 
1.17, and 1.4 THz in the refractive index and the absorption coefficient. 
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Figure 7. Optimum thickness for u.h.m.w. PE. The optimum thickness is determined from the measured absorption 
coefficient. The solid line is the second-order polynomial fitting to the absorption coefficient between 0.4 and 1.8 THz. 
The straight dotted lines indicate the thicknesses of 4.5, 5.5, 8, and 9 mm, available from the measurements. 
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Figure 8. Standard deviations in the optical constants of PE. Each standard deviation is determined from ten time- 
resolved signals. The scatter plots represent the raw profiles obtained from Equation 2, whereas the solid lines are from 
an analytical function, exp(aia; 5 + 02a; 4 + . . . + az,x + ae), fitting to the scatter plots. The result of the 9-mm pellet is 
very close to that of 8-mm, and thus is not shown here. 



the thicknesses of all the prepared pellets fall suboptimal at all frequencies of interest. Still the results can confirm 
that as the sample thickness approach optimum, the variance of the optical constants reduces proportionally. 

4.2.3. Standard deviation 

The standard deviations of the optical constants of PE at different thicknesses is shown in Figure 8. At all 
frequencies of interest, the standard deviation decreases, following the increment of the sample thickness. This 
is concurrent with the calculation of the optimum thickness, which shows the lowest standard deviation for a 
thicker sample. 

4.3. Liquid water 

The measurement for liquid water is different from that for solid dielectrics in terms of the propagation geometry. 
Water must be contained in a cell, equipped with a pair of transparent windows that maximise the throughput 
of T-rays. The present of the windows modifies the transfer function describing the measurement. Theoretically, 
the formula of the optimum thickness is still valid in this case. 

In the experiment, the cell in use is assembled from two parallel cycloolefines windows, each with the thickness 
of 3 mm, and a replaceable spacer in between, which permits the adjustment to the thickness of a liquid sample. 
Distilled water is injected into the cell, which is pre-adjusted to accommodate various thicknesses, from 15, 30, 
40, 60 (61), 80, to 170 //m. The sample with each thickness is measured with a focused THz beam at the normal 
angle of incidence for ten scans, with 45 sec separation between each. The reference is measured ten times with 
the identical cell, emptied with dry air, before and after the sample measurement. No significant change to the 
system and to the signal variance can be observed. 

4.3.1. Optical constants 

The optical constants of water, shown in Figure 9, are obtainable from the 170-/J,m sample measurement, which 
provides the lowest magnitude of Fabry-Perot fringes. Although normalised by the reference measured with 
a blank cell, the reflections cannot be completely removed. The phase at low frequency band is extrapolated 
from the phase value between 0.05 and 0.5 THz. The average refractive index is approximately 2.5, whereas 
the absorption coefficient rises linearly to 300 cm -1 at 1 THz. These values are in consistent with the optical 
constants of water at T-ray frequencies published in 12 and. 7 




Figure 9. Optical constants of water. The constants are calculated from the signal probing the 170-^im-thick sample. 
The available bandwidth of the measurement is from 0.05 to 0.8 THz. The apparent Fabry-Perot fringes are caused by 
reflections within the windows. 




Figure 10. Optimum thickness for water. The optimum thickness is determined from the measured absorption coefficient. 
The solid line is the second-order polynomial fitting to the absorption coefficient between 0.05 and 0.8 THz. The straight 
dotted lines indicate the thicknesses of 15, 30, 40, 60, 80, to 170 /im, available from the measurements. 
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Figure 11. Standard deviations in the optical constants of water. Each standard deviation is determined from ten 
time-resolved signals. The scatter plots represent the raw profiles obtained from Equation 2, whereas the solid lines are 
from an analytical function, exp(aix 5 + a2X A + . . . + a^x + ae), fitting to the scatter plots. 

4.3.2. Optimum thickness 

The exceptionally high absorption of water results in the requirement of a thin sample under measurement. 
Calculated from the absorption coefficient using Equation 10, the optimum thickness of water in the range 
between 0.05 to 2.0 THz is illustrated in Figure 10. It is clear that for this frequency range, the optimum thickness 
lies between 40 and 200 /im. Decreasing the thickness further than 40 /jm would cause the magnification of the 
variance. 

4.3.3. Standard deviation 

Figure 11 depicts the standard deviation of the refractive index of water measured at different thicknesses. 
Apparently, as the sample thickness decreases below 40 /im, the standard deviation increases proportionally at 
all frequency range of interest. At 1.7 to 2.0 THz, the 40-/mi sample shows the lowest standard deviation; at 
0.7-1.7 THz, the 60-^m sample; at 0.45-0.7 THz, the 80-/im sample; at 0.05-0.45 THz, the 170-//m sample. The 
results are in well agreement with the optimum thickness in Figure 10. The improvement of the precision can be 
clearly observed. By comparing the standard deviations of 15-/xm and 170-/xm samples at 0.2 THz, it is seen 
that the difference is nearly one order of magnitude. 

This experiment explicitly shows the impracticality of an excessively thin sample. Also it shows the function- 
ality of the proposed optimum thickness equation, although the transfer function of the sample in cell is different 
from that of the free-standing sample. 

5. USAGE OF THE MODEL 

Since the formula is frequency dependent, it can determine the thickness for any type of absorption response, 
provided that the absorption coefficient is available. In practice, this absorption coefficient can be measured 
from a preliminary sample or obtained from a published value. 

Two options are available, when the decision of the optimum thickness is made, either to have the widest 
measurement bandwidth or to have the lowest uncertainty at a particular frequency. If the widest bandwidth 
is required, the absorption value of material at the highest reliable frequency of the system should be used in 
determination of the optimum thickness. The sample with this thickness retains the bandwidth of the system, 
while provides reasonable uncertainty in measurement. For some reasons, one might need to observe the material's 
response at and nearby a particular frequency at the highest precision available from the instrument — for example 



a subtle absorption peak hidden beneath the noise. In this case, the optimum thickness calculated at that 
frequency of interest is no doubt the best selection. 

Despite the thickness optimality, for a very thin sample, the non-uniformities in the sample geometry and 
in bulk material tend to be pronounced. Combined with the effects from the frequency-dependent spatial 
distribution of the T-ray beam, relocation of the sample during measurements might cause a large variation 
among measured signals. 13 Although the variance in the signal's amplitude encompasses these effects, the 
dependence of the effects on the thickness is not taken into account in the optimisation. In addition, a thin 
sample magnifies the variance from thickness measurement. According to this concern, special care must be given 
to the uniformity of a sample with very thin optimum thickness and to the method of thickness measurement. 

6. CONCLUSION AND POTENTIAL EXTENSION 

Previously, in a transmission-geometry THz-TDS measurement, there is no explicit criterion for determining the 
appropriate thickness of the sample under measurement. The selection of the sample thickness depends wholly 
on the experience of an experimentalist. And in many cases, the sample thickness is preferred being as thin 
as possible, to achieve the widest bandwidth of the signal. This is favourable, unless one critical parameter is 
ignored — the uncertainty in measurement. An excessively thin sample can cause significant rise in the uncertainty, 
which is influenced by the noise in measured signals. 

This work offers a criterion in selecting the optimum thickness of a sample. Provided that the absorption 
of the sample at an interesting frequency can be estimated or approximated, the proposed model can predict 
the optimum thickness, which gives the lowest uncertainty in measurement. The derivation of the model is 
via minimising the uncertainty in optical constants — in terms of the variance or standard deviation — which is 
affected by the variance in measured time-resolved signals. 

The experiments, performed with three different materials, PVC, PE, and water, representing normal, low, 
and high absorption, confirm the validity and applicability of the proposed model. By selecting an optimum 
thickness for a sample, the improvement in the standard deviation of the optical constants as large as two order 
of magnitude can be observed. 

Validated by the experiments, the proposed model can be used as a rule of thumb in selecting the thickness for 
a sample to achieve the lowest measurement uncertainty. Furthermore, integration of this model into a frontend 
computer program, which provides on-spot measurement results, allows determining of an optimum thickness on 
site. 
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